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ABSTRACT 

We report the discovery of the Zeeman effect in the 44 GHz Class I methanol 
(CH3OH) maser line. The observations were carried out with 22 antennas of the 
EVLA toward a star forming region in OMC-2. Based on our adopted Zeeman 
splitting factor of z = 1.0 Hz niG^^, we detect a line of sight magnetic field of 
18.4 ±1.1 niG toward this source. Since such 44 GHz CH3OH masers arise from 
shocks in the outflows of star forming regions, we can relate our measurement 
of the post-shock magnetic field to field strengths indicated by species tracing 
pre-shock regions, and thus characterize the large scale magnetic field. Moreover, 
since Class I masers trace regions more remote from the star forming core than 
Class II masers, and possibly earlier phases, magnetic fields detected in 6.7 GHz 
Class II and 36 GHz and 44 GHz Class I methanol maser lines together offer 
the potential of providing a more complete picture of the magnetic field. This 
motivates further observations at high angular resolution to find the positional 
relationships between Class I and Class II masers, and masers at various frequen- 
cies within each category. In particular, CH3OH masers are widespread in high- 
as well as intermediate-mass star forming regions, and our discovery provides a 
new method of studying the magnetic field in such regions, by observing small 
physical scales that are not accessible by any other lines. 

Subject headings: ISM: clouds — ISM: magnetic fields — masers — polarization 
— radio lines: ISM — stars: formation 



1. INTRODUCTION 

The presence of magnetic fields in star forming regions presents theoretical and ob- 
servational challenges because of the computational complexity involved in incorporating 
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magnetic fields in simulations, and the difficult obs ervations that are required to measure 
their strengths and directions (ITroland et al.ll2008l ). While the role of magnetic fields i n 



regulating the onset of star formation is still a matter of debate (e.g., iCrutcher et al.ll2009l ). 
it has become increasingly clear that magnetic field s play a critical role in ca rrying angular 
momentum away from the protostar during collapse (IMcKee fc Ostrikerll2007l . and references 
therein). Moreover, the outflows along which this takes place may be driven by dynamically 
enhanced magnetic fields in the protostellar disk (IBanerjee fc Pudritzll2006l ). In particular, 
the driving of outflows along mag netic field lines may be crit ical in allowing accretion to con- 
tinue onto high mass protostars (IBanerjee fc Pudritzl 120071 ). Observational studies of such 
protostellar disks and outflows are difficult, due to the need for high angular resolution. In- 
terstellar masers, being compact and intense, offer a means of probing regions at high angular 
resolution. For years, their effectiveness as probes of star forming regions was overshadowed 
by a perception that the specialized conditions in which such masers form necessarily pre- 
vented them from being linked to conditions on larger scales. However, recent discoveries 
indicate that rather than being a measure in isolated atypical fragments , the magneti c fields 
measured in masers are indeed linked to the larger scale magnetic field. iFish fc ReidI (120061 ) 
found a relative consistency in the magnetic fields measured in clusters of mainline (1665 and 
1667 MHz) OH masers across a massive star forming region, and concluded that m a gnetic 
fields are ordered in massive star forming regions. More recently, IVlemmings et al.l (120 lOf ) 
have determined from polarization observations of 6.7 GHz CH3OH masers toward Cepheus 
HW2 that the masers probe the large scale magnetic field. 

The 6.7 GHz CH3OH masers mentioned above are examples of Class II methanol 
maser s, which are know n to be probes of the early phases of high mass star forming re- 
gions (jVlemmingsl l2008l . and references therein). Yet another class of masers, namely Class 



I me thanol masers, may probe even earlier phases of star forming regions (jPratap et al. 



20081 ). Such Class I CH. sOH masers are believed to be formed in collisional shocked regions 
in protostellar outflows (ICragg et al.lll992l : ISandell et al.l 120051 1. Together, Class I and Class 
II CH3OH masers offer a unique window into star forming regions by allowing us to observe 
the smallest physical scales that are not accessible via any other lines. Given the value of 
measuring the magnetic field in disks and outflows, we have embarked on a long-term plan 
to measure magnetic fields in CH3OH masers. In 2009, we discovered the Zeeman effect 
i n the 36 GHz Class I C H3OH maser line toward the high mass star forming region M8E 
( ISarma &: Momjianll2009l ). Th e Zeeman effect remains the most direct method for measuring 
magnetic field strengths (e.g.. ICrutcherl 119991 ). 



In this paper, we report the discovery of the Zeeman effect in the 44 GHz Class I 
CH3OH maser line. The observations were carried out with 22 antennas of the Expanded 
Very Large Array (EVLA) toward OMC-2. OMC-2 is a filamentary structure containing 
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several sites of active star formation (e.g., iTakahashi et al.ll2008[). and located about 12' 



northea st of the Trapezium OB cluster in Orion f lCastets fc Langerlll995l ). at a distance of 
450 pc fiGenzel fc Stutzkilll989l ). In §[21 we present details of the observations and reduction 
of the data. The analysis involved in extracting magnetic field information from the Zeeman 
effect is given in § [31 The results are presented and discussed in § [H 



2. OBSERVATIONS AND DATA REDUCTION 

Observations of the 7o — 6i A'^ methanol maser emission line at 44 GHz were carried 
out using 22 antennas in the D-configuration of the Expanded Very Large Array (EVLA) 
of the NRAOEI in two 2 hr sessions on 2009 Oct 25 and Nov 25. Table [1] lists the observing 
parameters and other relevant data for these observations. The data were correlated using 
the old VLA correlator, and in order to avoid the aliasing known to affect the lower 0.5 MHz 
of the bandwidth for EVLA data correlated with the old VLA correlator, the spectral line 
was centered in the second half of the 1.56 MHz wide band. The source 3C147 (J0542+4951) 
was used to set the absolute flux density scale, while the compact source J0607— 0834 was 
used as an amplitude calibrator. 

The editing, calibration, Fourier transformation, deconvolution, and processing of the 
data were carried out using the Astronomical Image Processing System (AIPS) of the NRAO. 
After applying the amplitude gain corrections of J0607— 0834 on the target source OMC-2, 
the spectral channel with the strongest maser emission signal was split, then self-calib rated 



in both phase and amplitude in a succession of iterative cycles (e.g.. lSarma et al.ll2002l ). The 
final phase and amplitude solutions were then applied to the full spectral-line uv data set, 
and Stokes I and V image cubes were made with a synthesized beamwidth of 1.93" x 1.58"; 
the maser emission was unresolved in these observations. Further processing of the data, 
including magnetic field estimates, was done using the MIRIAD software package. 



3. ANALYSIS 



For cases in which the Zeeman splitting Au^ is much less than the line width Ai>, the 
magnetic field can be obtai ned from the Stoke s V spectrum, which exhibits a scaled derivative 
of the Stokes / spectrum flHeiles et al.lll993l ). Here, consistent with AIPS conventions, / = 
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(RCP+LCP)/2, and V = (RCP-LCP)/2; RCP is right- and LCP is left-circular polarization 
incident on the antennas, where RCP has the standard radio definition of clockwise rotation 
of the electric vector when viewed along the direction of wave propagation. Since the observed 
V spectrum may also contain a scaled replica of the / spectrum itself, the Zeeman effect can 
be measured by fitting the Stokes V spectra in the least-squares sense to the equation 



V = al + - — 



b dl 
2 



(ITroland fc Heileslll982l : ISault et al.lll990[ ). The fit parameter a is usually the result of small 
calibration errors in RCP versus LCP, and is expected to be small. In these observations, 
a was of the order of 10~^ or less. While eq. ([1]) is strictly t rue only for thermal lines , 
numerical solutions of the equations of radiative transfer (e.g., iNedoluha &: WatsonI 119921 ) 
have shown that it gives reasonable values for the magnetic fields in masers also. In eq. ([1]), 
the fit parameter b = zB cos 9, where z is the Zeeman splitting factor (Hz n iG~^), B is the 



magn etic field, and 6 is the angle of the magnetic field to the line of sight ( ICrutcher et al. 



19931 ). The value of the Zeeman splitting factor z for CH3OH masers is very sma l l, bec ause 
CH3OH is a non-paramagnetic molecule. Following the treatment of I Vlemmings J2008h for 
the Zeeman splitting of 6.7 GHz methanol masers, we derive the Zeeman splitting factor 
using the Lande (^-factor based on laboratory measurements of 25 GHz methanol masers 
( iJenl Il95l[ ). and find z = 1.0 Hz mG~^. Since the actual Lande (7-factor for the 44 GHz 
line is not known, it is not possible to estimate the uncertainty in z. It is unlikely, however, 
that the measured Lande (^-factor will be vastly different from the value used here, since the 
-Bios values derived below seem to agree well with projections about the field and densities 
(see §11] below). In fact, if the measured value of the Lande (^-factor turns out to be vastly 
different, the Zeeman interpretation of these observations may have to be revised. Needless 
to say, these observations and the above considerations motivate immediate measurement of 
the Lande (7-factor in the 44 GHz CH3OH maser line. 



4. RESULTS AND DISCUSSION 

Figure [T] shows our Zeeman effect detection in the 44 GHz CH3OH maser toward OMC- 
2. As described in § [3l we determined magnetic fields by fitting the Stokes V spectra in 
the least-squares sense using equation ([1]). The fit parameter b (see eq. [1]) obtained by this 
procedure, together with the Zeeman splitting factor z = 1.0 Hz mG~^ discussed in §[3]above, 
then gives By^s = 18.4 ± 1.1 mC for the 44 GHz Class I CH3OH maser toward OMC-2. By 
convention, a positive value for i?ios indicates a field pointing away from the observer. 

In order to investigate the connection of our detected field to quantities in the larger 



- 5 - 



region, we looked at the 850 /im polarization map of iPoidevin et al.l (120 lOf ) observed with the 
JCMT (James Clerk Maxwell Telescope). Their map is reproduced in Figure[2l For reference, 
the position of our CH3OH maser is marked in this figure. We have used the position of 
the strongest 44 GHz CH3OH maser from the hi gh angular resoluti o n pha se referenced VLA 
A-configuration observations (HPBW 0^M5) of ISlvsh fc Kalenskiil (20091) . Associated with 



this position is a CO outflow (ITakahashi et al.ll2008l : IWilliams et al. 



20031 ) marked in red in 



Figure 121 the position angle of this CO outflow is 30°. Also show n in green in Figure [2] is 
the outflow detected in the 2.12 fim shocked H2 emission line from lStanke et al.l (120021 ): the 
position angle of this outflow is ~31°, identical to that of the CO outflow. In their high 



angular resolution VLA observations referred to above, ISlysh fc Kalenskiil ( 120091 ) observed 
six 44 GHz methanol masers spread out along a line aligned at an angle approximately 30° 
east of north. Upon compar ing all of the above information to the 850 /xm polarization 
map of iPoidevin et al.l (120 lOf ). it appears likely that the outflow is aligned along the local 
magnetic field, and the methanol masers are excited in collisional shocks along the outflow. 
Maser amplification is particularly efficient in directions approximately perpendicular to the 
shock propagation, and the compression of an ordered magnetic field from this orientation 
would give the following relationship: 



Bo 
Po 



Pi 



(2) 



where and 1 refer to the preshock a nd postshock (maser) regions respectively, and p is 



the gas density (e.g.. lSarma et al.l 120081 ). In the discussion below, we will use the molecular 



hydrogen number density n instead of p. If Bq and riQ are known, then we can use our 
detected value of i?ios to find n^ , a nd co mpare it to the density at which the 44 GHz CH3OH 
maser is excited. IPoidevin et al.l ( 120101 ) used the Chandrasekhar-Fermi method to estimate 
the value of the magnetic field in this region from their 850 pm observations, and found it 
to be 0.13 r nG: we adopt this as our value of Bq. We will use uq = 10^ cm~^, which was 



obtained bv ICastets fc Langerl (Il995[ ) based on C^^O observations toward OMC-2 (also see 
Poidevin et al.ll2010[ ). Finally, we adopt Bi = 2i?ios, based on statistical grounds (ICrutcher 



I999I ). Using these values in equation ([2]), we get rii ~ 10^ cm ^ for the postshock number 



density. This is in excellent agreement with theoretical models that s how the 44 GHz CH ^OH 
maser action is maximized in regions with density 10^ — 10^ cm~^ ( iPratap et al.ll2008l . and 
references therein). Another way of looking at this is that, knowing uq from observations 
and rii from theoretical models for CH3OH masers, our observed values of -Bios, and hence 
Bi, allow us to get a measure of the large-scale magnetic field Bq. 

Our detected value of i?ios in the 44 GHz Class I CH3OH maser line toward OMC-2 is 
similar to the fields (~ 25 mG) detected in our recent discovery of the Zee man effect in the 36 
GHz Class I CH3OH line toward the high mass star forming region M8E (ISarma fc Momjian 
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20091 ). It is also similar to t he fie lds obtained from the 6.7 GHz Class II CH3OH maser 



line observed by IVlemmingsl ( l2008l ). who detected significant magnetic fields with the 100 
m Effelsberg telescope in this line toward 17 sources, with an average value of 23 mG. The 
possibility of a fake Zeeman pattern due to beam squint is brought up every time there is a 
new Zeer nan detection, but i t is un likely in our maser observations; this is discussed in more 
detail in ISarma &: MomiianI (I2009I). Sin ce Class I and Class II CH3OH masers likely trace 
different spatial regions (jEllingserul2005l and references therein), the similarities in the fields 
measured in these masers suggests that CH3OH masers may trace the larger scale magnetic 
fields in star forming regions. Another possibility might be that Class I masers occur in the 
very early stages of star formation (before the formation of an ultracompact H II region), 
and Class II masers occur later on. In that case, the similarity in i?ios for these two classes 
may indicate that the magnetic field strength remains the same during the early stages of the 
star formation process. It is possible, however, that the similarities may result simply from 
selection effects due to orientation and/or the shock process. Therefore, a larger statistical 
sample as well as measurements of the Lande (^-factor are needed in order to test possible 
correlations or anti-correlations between fields measured in Class I and II masers and at 
different frequencies within each of these types. 

Finally, we use our observed value for the magnetic field to compare the magnetic 
and dynamical energies in these masing regions. The magnetic energy de nsity is given by 
3i?iQs/87r. Here, we have used the relation = 3Bf^^ that was obtained by ICrutcherl (119991 ) 
on statistical grounds. For i?ios ~ 18 mG from these observations, the magnetic energy 
density will be 4 x 10~^ ergs cm~'^. The kinetic energy density (thermal and turbulent) is 
given by {3/2)mna^, where a = Af /(81n2)^/^ is the velocity dispersion, the mass m = 2.8 nip 
(assuming 10% He; rrip is the proton mass). Since masers occur only in special directions 
along which they have developed the required velocity coherence, the velocity dispersions in 
the masing region may be great er than that t raced by masers. If we use Av = 5 km s~^ 



a reasonable value (see, e.g., lMcKeel[l999l ). and n = 10 cm , the density quoted in 



Pratap et al.l ( 120081 ) at which maser action in the 44 GHz line is maximized, we get a kinetic 



en ergy density equa l to 3 x 10 ergs cm ^. Even if Av is as high as 20 km s ^, as quoted 



by iTakahashi et al.l (|2008[ ) for the CO (3-2) outfiow marked at our CH3OH maser position 



in Figure HI we will get kinetic energy density equal to 5 x 10~® ergs cm~^. Therefore, the 
magnetic energy is dominant by at least an order of magnitude, indicating that the magnetic 
fields are dynamically significant in these regions. Also of interest is that the magnetic energy 
appea rs to be dominant in different regions supp orting Class I and Class II masers (e.g., this 
work; ISarma &: Momjiaru l2009l : IVlemmingj l2008l ) . 



- 7- 



5. CONCLUSIONS 

We have reported the discovery of the Zeeman effect in the 44 GHz Class I methanol 
(CH3OH) maser line. We detected a magnetic field -Bios = 18.4 ± 1.1 mC toward a star 
forming region in OMC-2. There may be systematic bias in this value due to the assumed 
Zeeman splitting factor of z = 1.0 Hz mG"^, which is based on laboratory measurements in 
the 25 GHz CH3OH maser line. Still, such values for i?ios imply that the magnetic energy 
may be dominant by at least an order of magnitude over the kinetic energy density, meaning 
that the magnetic field is dynamically significant in these regions. By assuming a feasible 
configuration for the shocks that give rise to 44 GHz CH3OH masers, we find that our 
detected value of -Bios leads to a value for the density in the postshock region (ni = 10^ 
cm~^) that agrees with the value predicted by theoretical models for maximizing maser 
action in the methanol line at 44 GHz. Alternatively, if we know the preshock density no 
(e.g., from molecular line observations) and ni (from theoretical models of masers), we can 
use our detected -Bios to find the large-scale magnetic fi eld. Together with the magnetic 



fields detected in the 36 GHz Class I CH3OH ma ser line f lSarma fc Momjiarul2009[ ). and 6.7 



GHz Class II CH3OH masers ( IVlemmingsl 120081 1 . our observations reveal the potential for 



obtaining a more complete picture of the magnetic field in star forming regions, since Class I 
and Class II masers are found at different distances from the protostar, likely trace different 
physical conditions, and may be excited at different stages of the star formation process. 

Given that Class I and II methanol masers allow observation of star forming regions on 
small physical scales that are inaccessible to other lines, and the possibility that these masers 
may trace the large-scale magnetic field, we believe that our discovery motivates several 
observational approaches. A high angular resolution study of the association between 36 
GHz and 44 GHz Class I methanol masers is now possible due to the new 36 GHz receivers 
on the EVLA. Likewise, a high angular resolution study of the association between Class 
I and Class II methanol masers is warranted. Identification of other maser sites where the 
Zeeman effect could be measured would allow for the conclusions above to be placed on a 
statist ical footing. Moreover , the OMC-2 region that we have studied is forming intermediate 



stars (ITakahashi et al.ll2008[ ). Observations of magnetic fields in 44 GHz CH3OH masers in 
high mass star forming regions, and the other studies described above, may yield important 
clues in the unsolved areas of high mass star formation. Finally, there is also a pressing need 
for laboratory measurements of the Zeeman splitting coefficients for various types of Class I 
and II CH3OH masers. 
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Table 1. Parameters for EVLA Observations 



Parameter 



Value 



Observation Dates 

Configuration 

R.A. of field center (J2000) 

Decl. of field center (J2000) 

Total Bandwidth 

No. of channels 

Channel Spacing 

Total Observing Time 

Rest Frequency 

Velocity at band center*^ 

Target source velocity 

FWHM of synthesized beam 

Line rms noise^ 



2009 Oct 25 & Nov 25 
D 

05^35°'27^66 
-05°09'39".6 
1.56 MHz 
256 
0.04 km s"^ 
4 hr 
44.069488 GHz 
13.2 km s"^ 
11.6 km s~^ 
1.93" X 1.58" 
P.A. = -10.40° 
8 mJy beam~^ 



'^The line was centered in the second half of the 1.56 
MHz band in order to avoid aliasing (see §[2]). 

'^The line rms noise was measured from the stokes / 
image cube using maser line free channels. 
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Fig. 1. — Stokes / {top— histogram) and V {bottom— histogram) profiles of tlie maser toward 
the 44 GHz Class I CH3OH maser in OMC-2. The curve superposed on V in the lower frame 
is the derivative of / scaled by a value of -Bios = 18.4 ± 1.1 mG. 



- 13 - 



I I I I I I I 1 T I I I 1 I T I I J I I r I T T r 




I ■ I I ) ' ' I I ■ * I I ' ' t I ' ' ■ I J I I 

-2 -4 -6 -8 -10 -12 -14 
La (drcmin,) 



Fig. 2. — Distribution of H2 jets (green) and CO outflows (red) in OMC-2 and QMC-3, super- 



i mpos ed on 850 fim polarization data (thin vectors), taken from Figure 7 of iPoidevin et al. 



(120101 ). The greyscale is the 850 /xm intensity, and is displayed to show the location of the 
filament. The reference position is 02000 = 05'^35™48^ ^2000 = — 05°00'0". The thick black 
cross shows the location of our 44 GHz CH3OH maser. 



